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Abstract: We present a detailed analysis of S-wave Kir scattering up to 2 GeV, making 
use of the resonance chiral Lagrangian predictions together with a suitable unitarisation 
method. Our approach incorporates known theoretical constraints at low and high energies. 
The present experimental status, with partly conflicting data from different experiments, 
is discussed. Our analysis allows to resolve some experimental ambiguities, but better data 
are needed in order to determine the cross-section in the higher-energy range. Our best 
fits are used to determine the masses and widths of the relevant scalar resonances in this 
energy region. 
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1 Introduction 



Chiral Perturbation Theory (%PT) [l]-0| provides a very powerful framework to study 
the low-energy dynamics of the lightest pseudoscalar octet. The QCD chiral symmetry 
constraints imply the pseudo-Goldstone nature of the it, K and r\ mesons and determine 
their interactions at very low energies. 

To describe the resonance region around 1 GeV, additional non-trivial dynamical in- 
formation is needed. One can construct a chiral-symmetric Effective Field Theory with 
resonance fields as explicit degrees of freedom g, [|. Although not as predictive as the 
standard chiral Lagrangian for the pseudo-Goldstone mesons, the resonance chiral La- 
grangian || turns out to provide interesting results, once some short-distance dynamical 
QCD constraints are taken into account ||. At tree level, this Lagrangian encodes the 
large- N c properties of QCD |10|,[ll|, in the approximation of keeping only the dominant 
lowest-mass resonance multiplets JT^ . 

The xPT loops incorporate the unitarity field theory constraints in a perturbative 
way, order by order in the momentum expansion. In the resonance region, the low-energy 
expansion breaks down and some kind of resummation of those chiral loops is required to 



satisfy unitarity [pT^ -|2T|. This effect appears to be crucial for a correct understanding of 
the scalar sector, because the S-wave rescattering of two pseudoscalars is very strong. 

In this paper, we present an attempt to describe S-wave Kir scattering up to 2 GeV, 
making use of the resonance chiral Lagrangian predictions together with a suitable unitari- 
sation method. The analysis is far from trivial, owing to the very unsatisfactory experi- 
mental status of the scalar sector, with partly conflicting data from different experiments. 
Even the masses and widths of the dominant scalar resonances are not clearly established, 
since the large unitarity corrections tend to mask the resonance behaviour. 

A good understanding of the scalar isospin-1/2 system is required in order to improve 
the determination of the strange quark mass from QCD sum rules for the strangeness- 



changing scalar current [ 22|— p7] , since the phenomenological parametrisation of the corre- 
sponding spectral function can be related to the Kn scalar form factor. We plan to discuss 



this form factor in detail in a subsequent publication [28 



In section 2 we collect the one-loop %PT results for the Kn scattering amplitude, 
including scalar and vector resonance contributions. Our unitarisation approach for the 
elastic channel is discussed in section 3. Section 4 analyses the experimental data up to 
1.3 GeV and several xPT-inspired fits are performed. At higher energies, one needs to 
include inelastic channels; the dominant two-body Ktx — > Krj, Krj' modes are incorporated 
in section 5, and used in section 6 to fit the higher-energy data up to 2 GeV. Our fits 
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determine the resonance poles of the T matrix in the complex plane; this is worked out 
in section 7, where we compare results for the resonance masses and widths from different 
fits. Our conclusions are briefly summarised in section 8. To simplify the presentation, we 
refer to refs. [f§-|9| for details on the chiral formalism and notations, while some technical 
aspects and cumbersome formulae are relegated into appendices. 



2 Chiral expansion for the Kn scattering amplitudes 

Up to one loop in the chiral expansion the Ktt scattering amplitudes have been calculated 
by Bernard, Kaiser and MeiBner |29| and resonance contributions have been included by 



the same authors in fgOj. We shall just review the main expressions for the Ktt amplitudes 
to make our work selfcontained. 

Since the pion and the kaon have isospin 1 and 1/2 respectively, there exist two inde- 
pendent Kit amplitudes T\l 2 and T%/£ (in the following, for simplicity, we shall drop the 
subscript "Kir"). The charge-two process K + tt + — > K + ir + is purely 1 = 3/2 whereas the 
process K + tt~ — ► K + tt~ contains both 1 = 1/2 and 1 = 3/2 components. Both amplitudes 
depend on the Mandelstam variables s, t and u and are related by s <-> u crossing. From 
this relation one finds[] 

T 1/2 (s, = | T 3 / 2 (n, t, s) - \ T 3 /\s, t, u) . (2.1) 

Therefore it is sufficient to give here only the 1 = 3/2 amplitude T 3 / 2 . 

At order p 4 in the xPT expansion including resonances, T 3 / 2 can be decomposed as 
follows 

T 3/2 = T (2) + T (4) + T (4) + T {4) ^ (2 _ 2) 

where is the leading order p 2 contribution already known from current algebra |32] . 



l£\ T< 4 > and If 

are order p A resonance, unitarity and tadpole corrections respectively. 
We have assumed that at the resonance-mass scale the low-energy constants in the order 
p 4 chiral Lagrangian are saturated by resonance exchange @,[£j. The tree-level contribution 
T ( 2 ) is given by 

T {2) (s) = -}—(M* K + Ml-s). (2.3) 

^JKJtt 

In writing eq. ( |2.3| ), we have absorbed some order p 4 corrections into by replacing l// 2 

with 1/UkU)- 

1 For further details the reader is referred to refs. [EMKlj]. 
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Below, we list explicit expressions for the order p 4 contributions , and T^ 1 to 
the amplitude T 3//2 in %PT with resonances. The resonance contribution T^\s,t,u) is 
given by 

) u(t - .< 



T, 



(i) 



r<2 



R 



2 f2 

J IT 



IK 



t[U — S) 



A 2 



u 
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c d u 



Ml ^ 



1 

+ 3 



M Si 



Mg — t 



c d t-2{c d -c m )M 2 K c d t-2{c d -c m )Ml ,(2.4) 



where we defined the constants A Kw = {M\ — M 2 ) and J2 Kn = [M\ + M 2 ). Gy as well as 
c d and c m are coupling constants of the vector and scalar mesons in the chiral resonance 
Lagrangian of refs. ||[| respectively. For the coupling constants of the scalar singlet meson 
Si, c d and c m , we have already used the large- N c constraints c d = c d / \/3 and c m = c m /\^3. 
Ms refers to the mass of a generic scalar octet meson and we shall always employ the 
large- N c result M Sl = M s . 

The unitarity correction Tjj\s, t, u) has originally been calculated in ref. |29[] and takes 
the form 

1 
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tiu — s , 
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(2.5) 



We have already corrected for two misprints which appeared in the original publication |29 



and the explicit form of the loop functions can be found in appendix A. The unitarity 
correction has also recently been checked independently in ref. ||33|| . 
Finally, the tadpole contribution reads 

A K tt 



Tr 



(1) 



3«tt ~ — «, 



with fip being functions which depend on the renormalisation scale u and are given by 



2.6) 



M - Inl 



Mf 



(2.7) 
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It has been assumed that the resonance contribution saturates the low-energy constants 
Li and thus, in accordance with ref. ||, in the following we always set the renormalisation 
scale \x = M p , where this saturation is found to take place. Both, the resonance as well 
as the tadpole contribution above, take different forms compared to the original papers 
|29|,|30[], because we chose to write the tree-level term ( |2.3|) differently However, we have 
independently calculated these terms finding agreement with [[SOf . 

Since in this work we are only interested in the S-wave scattering amplitude, the 
partial wave projection still has to be performed. Generally, the partial wave amplitudes 
t{ (s) are given by 

1 ) 

ti(s) = — JdzP^T^z), (2.8) 

-l 

where / is the total angular momentum, z = cos 8 is the scattering angle in the cm. system 
and Pi(z) are Legendre polynomials . 



3 Unitarisation of the amplitudes 

A numerical inspection of eq. ( |2.2| ) shows that already at energies of the order of 1 GeV the 
unitarity correction is rather large. Since we aim at describing S-wave Kir scattering 
well above 1 GeV, we have to resort to methods which unitarise the partial wave ampli- 
tudes. The approach we are going to use is based on the N/D method |35| and has been 
developed in refs. [p0|,pl],[36[|. It can be derived by working out the partial wave amplitudes 
from a perturbative loop expansion in terms of the unphysical cuts. The generated partial 
waves are then matched with the ones derived from next-to-leading order xPT @ including 
in addition the explicit exchange of resonances |§ . 

Restricting oneself to two-body unitarity and neglecting the unphysical cut contribu- 
tions, with the help of the N/D method [35] in ref. 0] the corresponding general structure 



of a partial wave amplitude was derived. This structure was shown to match with the 
local tree-level %PT amplitudes plus the s-channel exchange of explicit resonance fields. 
It was also estimated that the contributions of the unphysical cuts in the physical region 
below 1 GeV for the 1 = 1/2 meson-meson S- waves amount just to a few percent. Then, 
in ref. |21],[36| the generalisation of this approach to include the unphysical cuts up to one 
loop in x?T was established. Note that the present expansion, treating the unphysical cut 
contributions as small quantities, can be implemented as a chiral loop expansion of these 
contributions. This means that at the tree level our amplitudes are crossing symmetric, 
e.g. resonance exchanges are included in all the s, t and w-channels, and that the resulting 
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partial waves amplitudes fulfil the N/D equations exactly in the right hand cut and up to 
one loop in %PT for the unphysical cuts. This is precisely the order required to match 
with the chiral expressions of ref. where the local and pole terms were supplied 
with loops calculated at 0(p 4 ) in xPT. In refs. @, 0] the method described above has 
already been used to study the strongly interacting Higgs sector and elastic tcN scattering, 
respectively. 

Up to the considered order, we can write the unitarised partial wave amplitude t\ (s) 
as 0:0 

tj {s) = — *M — (3.i) 

with 

N{{s) = t I l (s)-g I Kn (s){t I l {2 \s))\ (3.2) 

and 

9Kn( s ) = 16vr J^(s) + c^ ff . (3.3) 

On the right-hand side of eq. ( |3.2|) , (s) refers to the leading 0(p 2 ) expression of tj(s), 
with tj (s) itself being calculated at the one-loop level in x?T plus resonances. J(s) is the 
standard two-particle one-loop integral as given in appendix A, and the c I Kn are arbitrary 
constants, not fixed by the unitarisation procedure, which have to be real in order to satisfy 
unitarity. Reexpanding ( |3.1| ) up to one loop, also taking into account eq. fl3.2|), one can 
immediately match the resulting expression to the chiral amplitudes collected in section 2. 

It is easy to see that eq. (fOJ satisfies elastic unitarity above the Kn threshold and 
below the inelastic ones, namely: 

Imtj(s) = a X7r ( S )|t7( S )| 2 , (3.4) 

since by construction ImiV/(s) = in the Kir physical region and 



Img^s + iO) = a Kn (s) = yj(l-s + /s)(l-s_/s), (3.5) 

for s > s + where s + = (M K + M w ) 2 and s_ = (M K - M^) 2 . 

As long as we only consider a single channel, tj(s) and Nf(s) are scalar functions. In 
ref. |n[] the rather straightforward generalisation of eq. ( |3.1| ) to coupled channels where 
tj(s) and Nf(s) become matrices is also given and we shall come back to it in section 5 
below. 



2 For more details, we refer the reader to section 3.2 of the review [Ell 
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4 Fitting the elastic channel 



Experimentally, it has been observed that the 1 = 1/2 Kir scattering amplitude turns out 
to be elastic below roughly 1.3 GeV and the same was found for the I = 3/2 amplitude up 
to 1.9 GeV. Therefore, in this energy range we should be able to describe the experimental 
data on Kir scattering with the unitarised chiral expressions presented in the previous two 
sections. First we shall comment on the experimental data which we included in our fits, 
before we explain the actual fitting procedure and our results. 

The two most recent experiments on Kn scattering have been performed by Estabrooks 
et al. |38| and Aston et al. [|3~9"1 , fjQ|. Estabrooks et al. measured charged Kir scattering 
for all four charge combinations and were thus able to separate the I = 1/2 and 7 = 3/2 
components in the elastic region below 1.3 GeV. Aston et al. only considered K~tt + 
scattering. In this case the S-wave amplitude in terms of phase shifts and inelasticities is 
given as follows: 

A = a c *> = AT + \ AT = \ % (vl^ - l) + ± (r&^ - l) , (4.1) 

where t]q and 5q are the inelasticities and phase shifts for the 1 = 1/2 and 1 = 3/2 channel 
respectively. In writing eq. (|4.1|) we have adopted the normalisation of refs. [|38| -|l0|| f° r A)- 
In ref. |j38| , 5^ 2 has been extracted for = 0.73 - 1.30 GeV, and S^ 2 for 0.73 - 1.72 GeV. 
In this work it was also demonstrated that below 1.3 GeV, the scattering is purely elastic, 
that is, r^ 2 = ifj 2 = 1, whereas it was found that rfj 2 = 1 in the full energy range. On 
the other hand in ref. |4(J only a$ and 0o up to \/s = 2.52 GeV were given. 

Because the 1 = 3/2 component of the scattering amplitude was only measured directly 
by Estabrooks et al. [[38], we have decided to also include older data for 8q in our fits 
J|lH44l. The last three experiments only presented data for the K~tt~ cross section and one 



Q /2 

still has to extract 5 from the given cross sections. In addition, rather different energy 
intervals have been used by the experimental groups |JS|,[IT]-0 for calculating cross sections 
and phase shifts. To be able to combine the data consistently, we thus also included an 
error in the energy. Nevertheless, it turned out that the combined x 2 /d-o.f. for our fits 
was only of the order of one if we increased the error in all the experiments [p^,|4T|-^4| by 
a factor of two. The increase of errors for these data is implied in the rest of our work. 
In terms of the unitarised chiral expressions, the isospin amplitudes Aq are given by 

K(s) = ? ^MMM (4.2) 

In the elastic region the phase shift 5q is then related to the amplitude by tan^Q = 
Iul4q/IW1q. We have fitted all data for 5^ 2 , a and 0o up to 1.31 GeV and for 5q 2 
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up to 1.9 GeV simultaneously. The parameters for the fit are Mk*, cj, c m , cj^ and c%^. 
Lacking further information, and to reduce the number of free parameters, we have set 
M Sl = M s = M K *, and we have used the large- N c estimate G v = = 65.3 MeV [0. 

The latter value is in good agreement with the lowest-order results coming from the decay 
width of the K* (p) meson which yield Gy = 65 (67) MeV and is also compatible with the 
original estimate Gy = 53 MeV || which however already includes chiral loops. All other 
parameters have been taken from the Review of Particle Physics [[45] . For the convenience 



of the reader we have compiled the constant input parameters in appendix D. 

From our fits we observe that some of the parameters are strongly correlated and for 
example rather different values for the couplings of the scalar resonances can lead to similar 
fit quality. Thus it is necessary to have an idea what the values of Cd and c m should be. 
In ref. H values for these parameters were estimated by assuming resonance saturation 
for the low-energy constants in the chiral Lagrangian. Taking the same approach with the 
most recent values for the low-energy constants [46, [471, and an average mass of the scalar 
octet Ms = 1.4 GeV, we arrive at 

\c d \ = 30 ±10 MeV, \c m \ = 43 ±14 MeV, (4.3) 

with Cd c m > 0, rather close to the original estimate in ref. ||. For the error we have taken 
one third which is a typical error for large- N c estimates of chiral constants. The estimate 
Ms = 1.4 GeV tacitly assumes that the correct 1 = 1 scalar in this region is the ao(1450), 



and that the a (980) is generated dynamically [|lq , |20|, |4q-pTl ■ Note however, that this 
question is controversial in the literature |5*^|-p^] ,p| 

As a first fit, we fix q and c m to the central values presented in eq. ( |4.3j ) above. The 
remaining fit parameters are then found to be 

M K * = 1.29 GeV, c^ = -0.094, cf n = 1.703, (4.4) 

with a x 2 — 223.5/117 d.o.f. = 1.91. The fits are performed with the program Minuit p5[ . 
Since the fit parameters are highly correlated, we do not give statistical errors for the fit 
because they would underestimate the true uncertainties. We shall return to a discussion 
of the uncertainties below. 

Next, the Kir scattering data are fit with the constraint c m = Cd for the couplings of 
the scalar resonances. This constraint is motivated from a calculation of the scalar form 
factor including resonances in the large- iV c limit which will be presented in a subsequent 



3 For further discussion and references on this point see also the note on scalar mesons in the Review 
of Particle Physics M5[ . 
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publication ||28|| . Requiring the tree- level form factor to vanish at infinity, which should be 
a plausible assumption, in the case of one scalar resonance we find that q = c m = f n /2. 
For the corresponding fit the parameters take the values 



M K * = 1.19 GeV, 



Cd 



c 1/2 



= 45.4 MeV. 
1.689. 



-0.376, c%* = 1.689. (4.5) 

The x 2 f° r this fit is 210.6/116 d.o.f. = 1.82. The resulting values of q and c m turn out to 
be surprisingly close to the large- N c constraint f w /2 = 46.2 MeV. 

Finally, we perform a fit leaving both q and c m as free parameters. This fit results in 



M, 



h~* 



1.26 GeV. 



.1/2 



Cd 



24.8 MeV. 



76.7 MeV 



0.191 



.3/2 



1.700 



(4.6) 



and x 2 — 201.1/115 d.o.f. = 1.75. Here, c m turns out to be quite different from the central 
estimate given in eq. ( |4.3|) . However, comparing eqs. ( |4.4|) , (|4.5| ) and (|4.6| ) one sees that 
the x 2 is around two for all three fits. Thus, the difference for the estimates of q and c m 
gives an indication about the systematic uncertainties in their determination. 

The fit result corresponding to eq. fl4.6|) is shown as the long-dashed line in figs. 1-4 



for 5, 



1/2 



ao, 



and 5 respectively as a function of the invariant mass of the Kir system 
Mku — y/s- The curves corresponding to the other two fits look similar. In figs. 2 and 3, 
data for ao and 0o up to 1.6 GeV has been included already. The notation for the different 
data sets is as follows: ref. R8J full circles; ref. [03,^01] full squares; ref. HTJ] full triangles; 



ref. E| open circles; ref. E3J open squares; ref. MM open triangles. As discussed above the 



errors for the data by refs. [38,41-4_4J have been increased by a factor of two. Concerning 
the quality of our fits we observe that the total x 2 1 d.o.f. is around two. This is due to 
the fact that the x 2 f° r ^o'' 2 individually is almost three. The individual x 2 °f ^c/ 2 \ a o an d 

3 /2 

</>o are all close to one. Even removing the data for % above 1.3 GeV only reduces the 
X 2 for 5q^ 2 to about two. Thus we conclude that the unitarised chiral expression does not 
provide a fully satisfactory fit to the 5^ 2 data and that it somewhat prefers the older data. 
Nevertheless, compared to the conventional xPT result without unitarisation p9[[3(J, the 
situation is much improved. For a comparison with the inverse amplitude method (IAM) 



see also ref. 56 



The somewhat unsatisfactory description of the 1 = 3/2 data is probably due to the 
fact that there is no resonance contribution in the s-channel. Being pure background the 
chiral expressions are not expected to work very well at higher energies, because additional 
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Figure 1: Elastic fits for o . Full circles: data of ref. ||38|| . Long-dashed line: fit of 
eq. (|4.6| ); solid line: fit of eq. ( |4.10| ); short-dashed line: fit of eq. (|4.13| ); dotted line: fit of 
eq. O. 

contributions not included in our low-energy description (e.g. higher-mass resonance ex- 
changes in the t and u channels) could be important numerically. This is reflected in the 
fact that the contribution from unphysical cuts is much larger in the case of 7 = 3/2 than 
for 7 = 1/2. However, the discrepancies among different sets of data, and the fact that the 
older 1 = 3/2 data admits a better theoretical description, prevents us from drawing any 
clear conclusion. More precise experimental data are needed in this channel. 

To further improve our overall fit for all data sets, we decided to use a unitarised 
background parametrisation for the 5q^ 2 data: 



tan Sq^ 2 (s) 



aq Kv {s) 



(i + /*gL(«)+7?Us: 



(4.7) 



where qKn(s) = y/s <j j R- 7r (s)/2 is the Kn centre-of-mass momentum and a, (3 and 7 are 
parameters of the ansatz. The parameter a can be identified with the scattering length. 
Again, we perform the same three fits as with the chiral expression for the 1 = 3/2 channel. 
First fixing a and c m to the central values of eq. Q4.3|) , we obtain 
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Figure 2: Elastic fits for a^. Full circles: data of ref. []38|1 ; full squares: data of ref. l39],|40" 
Long-dashed line: fit of eq. ( |4.(j| ); solid line: fit of eq. ( |4.1U| ); short-dashed line: fit of 
eq. (P^D ; dotted line: fit of eq. flOg) . 



« = -0.897 GeV" 1 , (3 = -0.284 GeV" 2 , 7 = 1.945 GeV" 4 , (4.8) 

with a x 2 = 107.7/115 d.o.f. = 0.94. The fit parameters imposing the constraint c m = Cd 
are found to be 

M K * = 1.19 GeV, Cd = c m = 44.8 MeV, cj? = -0.381, 

a = -0.859 GeV -1 , (3 = -0.580 GeV" 2 , 7 = 2.301 GeV -4 , (4.9) 
and x 2 = 96.4/114 d.o.f. = 0.85. Finally, leaving both Cd and c m as free parameters yields 

M K * = 1.36 GeV, c d = 13.0 MeV, c m = 85.4 MeV , c^J = 0.200, 

a = -0.860 GeV" 1 , (3 = -0.567 GeV" 2 , 7 = 2.283 GeV" 4 , (4.10) 

with x 2 = 90.5/113 d.o.f. = 0.80. The last fit is shown as the solid curve in figs. 1-4. 
Compared to the fits with the chiral expression for the 1 = 3/2 amplitude, the total 
X 2 / d.o.f. is greatly improved. The individual x 2 f° r ^c/ 2 i s around one and all other x 2 
are below one. The parameters a, (3 and 7 have now replaced the parameter c 3 ^ 2 in the 
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Figure 3: Elastic fits for cj) . Full circles: data of ref. p8| ; full squares: data of ref. p9|.[40 
Long-dashed line: fit of eq. (|4.6| ); solid line: fit of eq. (|4.10|) ; short-dashed line: fit of 
eq. flOSl ); dotted line: fit of eq. 



original parametrisation of the 1 = 3/2 amplitude. We observe that for the fits with q 
unequal c m , the Kq mass turns out to be somewhat larger than the results in eq. ( |4.5| ) and 
( (4.9| ) and the value of depends very much on the other parameters. However, again 
rather different values of c<j and c m lead to similar \ 2 implying sizeable uncertainties for 
these parameters. 

Next, we compare our fits for the scattering amplitudes with the unitarised chiral 
expressions to a parametrisation for the I = 1/2 amplitude with a .fT-matrix ansatz pi , 
57, keeping the background parametrisation of eq. ( [4.7| ) together with the parameters 
of eq. ( 4.10 ) for the 1 = 3/2 amplitude. The conventional f^-matrix ansatz is completely 
analogous to eq. (|4 



a K Js)K, 



1/2, 
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1 /2 

where now C(s) = 16ttJktt(s), that is 



- C(s) K 
0, and 



1/2, 
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K, 



1/2, 



+ 
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{Ml-s) ' (1 + cs) ■ { ■ ] 

The i^-matrix ansatz consists of a resonance term, with gp> being the coupling of the 
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Figure 4: Elastic fits for 8q . The notation for the data sets is as follows: ref. p8| full 



circles; ref. [4l[ full triangles; ref. |42| open circles; ref. |43[ open squares; ref. J34| open 
triangles. Long-dashed line: fit of eq. flOD ; solid line: fit of eq. (|4.10| ); dotted line: fit of 
eq. flQ4]) . 



resonance, and a background contribution. We have tried different forms of the background 
term, but the fit does not depend very much on this choice. In figures 1-3 the i^-matrix fit 
corresponds to the short-dashed lines and the corresponding fit parameters take the values 



M 



R 



1.32 GeV. 



0.695 GeV" 



9r = 0.628 



c = 0.172 GeV" 



(4.13) 



The x 2 = 86.9/116 d.o.f. = 0.75 for this fit. Thus the fit quality is very similar to our 
fit (|4 . X Q|) with the chiral expressions where we also used a pure background form for the 
1 = 3/2 amplitude. 

Nevertheless, as can already be anticipated from the curvature of the two fits close 
to the Kit threshold, the I = 1/2 scattering lengths turn out to be somewhat different. 
Expanding our results close to threshold, we obtain ay 2 (K-matrix) = 0.31 M" 1 , whereas 



a 



1/2 



(chiral) = 0.18 M^ 1 , although the quality of the fit for / = 1/2 is equally good. 
As has become conventional in the literature, we have expressed the scattering lengths 
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in units of M~ x . Even more dramatically, for the 1 = 3/2 scattering lengths we find 
a^ 2 (background) = — 0.12M" 1 and ag^chiral) = — 0.05 M" 1 . Our %PT results for the 
scattering lengths are in agreement with the order p 4 calculation of ref. [P| (table 3). 
These findings show that reasonable fits to the experimental data can be obtained that 
spoil the xPT behaviour in the low-energy region close to the Kn threshold. We conclude 
that the scattering lengths for S-wave scattering can only be extracted more reliably 
once better experimental data close to threshold become available. 

Finally, to get some feeling about the importance of the resonance contribution and the 
real part of the unitarity corrections, as the dotted line in figs. 1-4, we display the leading 
order x;PT result using a unitarisation according to the prescription: 

tan^ = <x AV ( S )^ (2) ( S ). (4.14) 

This prescription is equivalent to eq. (|3.1|) with the choice g(s) = iaKn(s) and Nq(s) = 
t (s). We should remark that this result is a parameter free prediction of lowest-order 
%PT. However, figs. 1-4 show that, except for 0o below roughly 1.2 GeV, resonance and the 
real part of unitarity corrections are required in order to obtain an acceptable description 
of the data. 



5 Including the Krj and Krj f channels 

Above roughly 1.3 GeV, Kn scattering in the 1=1/2 channel was found to be inelastic 
[38|,[39[] with Krj' being the first important inelastic channel. Therefore, in order to be able 
to describe Kn scattering beyond 1.3 GeV, we need to include the Krj' channel in the 
scattering amplitude in the framework of xPT with resonances. Although the Krj channel 
was found to be unimportant, nevertheless, for consistency we have also included it in our 
analysis. The importance of the Krj channel for our fits shall be further discussed below. 

In the large- A^ c limit, the singlet rji field becomes the ninth Goldstone boson and can 
then be incorporated with an extended U(3)l®U(3)r chiral Lagrangian The large 

mass of the rf originates from the U(1)a anomaly, which, although formally of 0(1/N C ), 
is numerically important and cannot be treated as a small perturbation. It is possible to 
build a consistent combined chiral expansion in powers of momenta, quark masses and 
l/N c , by counting the relative magnitude of these parameters as m q ~ l/N c ~ Oijp 2 ) |^TJ. 

The physical states rj and rj' are mixtures of the SU(3) singlet and octet states rji and 
7/g . If isospin-breaking effects are neglected, the mixing can be parametrised in terms of 
the mixing angle 9 usually defined by 
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1 77) = cos 6* \r]s) — sin 9 \rji) , 
\r]') = sin 9 \r] 8 ) + cos 9 \rji) . 



(5.1) 



For our numerical analysis below, we shall use sin^ = —1/3 ~ —20° which is in the 
ball-park of present day phenomenological values. Further details on 77-7/ mixing in the 
framework of %PT with resonances can be found in appendix B. To be able to include 
the Kr) and Krj' channels in the analysis of Kir scattering we have to calculate the five 
scattering amplitudes Tk-k^k^ and Tk^^k^ with i,j = 1,8. The resulting expressions 
are rather lengthy and therefore, although representing one of our main results, have been 
relegated to appendix C. Eqs. ( C.l ) to ( |C.5| ) include the leading order p 2 and the resonance 
contributions. So far, we have not calculated the order p 4 loop contributions to the inelastic 
channels, but as will be explained in more detail below, the dominant contributions are 
accounted for by our unitarisation procedure. Using eq. Qj.ID , the five required amplitudes 
Tktt^Kti, TxTr-tKri', Tk^ktji Tk^Kt) 1 an d Tk^'^Kt, 1 can then easily be expressed in terms 
of T K7T ^ Krj . and T Klli ^ Krj] . 

For a coupled channel analysis, our unitarisation procedure of section 3 has to be 
generalised and in particular eq. ( |3.1| ) has to be modified to a matrix equation. In the 
three-channel case, N(s) and g(s) become three by three matrices. Since in the following 
only S-wave Ktt scattering in the 1 = 1/2 channel is considered, to simplify the notation 
we have dropped the corresponding indices in the expressions. Then the generalisation of 
eq. fl3?Tj) takes the form [|~ 
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t(s) 



N(s) 



g(s)N(s) 



with 



N(s) 



( N K k-+k-k{s) N K7r ^ Kv (s) N K7T ^ Kv/ (s)\ 



N. 



Kt]^Kri\S) 



N. 



Kri->Kri' 



(5.2) 



(5.3) 



and 



\N K ^ Kr f(s) N Krj ^ Kv ,{s) N Kr)l _ Kri f(s) J 



(5.4) 



/WOO 

g(s) = 9k v {s) 

V g Kv '{s)J 

Several remarks are in order at this point. Because of time-reversal invariance, N(s) has to 
be symmetric. In addition, to satisfy unitarity (s) needs to be real above the highest 
threshold of the channels i and j, and 



Im gi 



(5.5) 
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where i denotes the three channels Kn, Krj or Kr]', Sj is the corresponding threshold, and 
the cTj(s) are defined in analogy to eq. ( |3.5|) . 

Nk^Kk is already given by eq. ( gg) , and N KlT ^ Kri , N Kn -> Kv >, N K ^ Kv , N Krt ^ Krjl as 
well as NKr) r -*Kr]' result from the S-wave projection of the corresponding T amplitudes 
according to eq. ( p.8|) . No subtractions have to be performed for these N(s) because we 
have not included the loop contributions to the Kr\ and Krj' channels explicitly. Therefore 
in this case N(s) is actually equal to t(s) (compare with eq. ( |3.2|) ). Finally, analogous to 
eq. (|3.3|) , gKn( s ) an d gi<r]>(s) are written as 

g K r,(s) = 16irJ Kv (s) + c Kr] , 

gKrf{s) = 16iiJ Kv ,(s) + c Kv < , (5.6) 

which satisfies eq. ( |5.5|) automatically. Again, ck v and c^y are arbitrary real constants. 

As it stands, eq. (|5.2|) already includes the dominant part of the unitarity corrections 
which arise from the loop diagrams in x?T also for the inelastic channel. The imaginary 
parts are generated by the imaginary parts of gK V (s) and gK V '(s). We have verified ex- 
plicitly that the imaginary part generated from the leading contribution to the T^-^k^s 
amplitude agrees with the imaginary part of the Tktt^k-k resulting from the r] 8 terms of 
eq. ( |2.5|) . Thus, when using eq. (|2.5| ) in our numerical analysis, for consistency we have 
dropped the rj 8 contribution since it arises from the Tk^k^ amplitude. On the other 
hand, the real parts of the unitarity corrections in the inelastic channels are not complete; 
the missing pieces are parametrised by the constants ck v and c^y. 



6 Fitting the inelastic case 



For the inelastic case, the unitarised chiral expression corresponding to the S-wave I = 1/2 

1/2 

amplitude A takes the form 

-l 



<'\s) 



a(s)N(s) 



g(s)N(s) 



cr (s 



(6.1) 



where 

'vkM \ 

a(s) = a Kv (s) . (6.2) 

V (T Kv ,{s)J 

The Kit — > Kit contribution is obtained by taking the relevant matrix element of the 

1 /2 

matrix A Q (s). If the Kr\ and Krj 1 contributions are removed one of course recovers the 
corresponding expression (|4.2|) for the single channel case. 
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Let us next discuss the situation of the experimental data in the inelastic region above 
1.3 GeV. Because experimentally only the absolute value of the scattering amplitude can 
be determined, there are 2 L discrete ambiguous solutions to the partial wave decomposition 
of the data where L is the maximal angular momentum that contributes at a given energy. 
These discrete ambiguities are apparent from the behaviour of the imaginary part of the 
amplitude (Barrelet) zeros Zi |62], if the scattering amplitude is expressed as a complex 
polynomial in z = cos 8 where 9 is the scattering angle: 

A(s,z) = X)(2Z + l)A,( a )P I (^) = f(s)f[(z-z i (s)). (6.3) 

1=0 i=l 

The measurement of |A(s,z)| cannot determine whether Zi(s) or z*(s) is a zero of A(s,z); 
thus the 2 L ambiguity arises. Requiring the solutions to be smooth, it is possible to switch 
from one solution to another when a Barrelet zero Z{ approaches the real axis. In the 
elastic region the ambiguity can be fixed from unitarity. However, in the inelastic region 
the behaviour of the Barrelet zeros has to be studied in detail. 



Whereas in the experiment by Estabrooks et al. 0] a fourfold ambiguity arose above 
1.5 GeV, because the imaginary parts of the two Barrelet zeros z\ and Z2 were found to 
vanish within the experimental uncertainties, the group by Aston et al. |3S|] obtained an 
unambiguous solution up to 1.85 GeV. In the region between 1.86 GeV and 2.0 GeV 
Aston et al. found two solutions whereas above 2.0 GeV four solutions to the partial wave 
decomposition remained. Inspecting the signs of the imaginary parts of the z^ it is clear 
that the solution of Aston et al. below 1.85 GeV corresponds to solution B of Estabrooks 
et al. Since both experiments have similar statistics in our opinion the question whether 
the solution of Aston et al. in the region between 1.5 GeV and 1.85 GeV is the physical 
one, requires further corroboration. Nevertheless, we have decided to perform our fits to 
the experimental data by Aston et al. |39|,^(| and solution B of Estabrooks et al. [38] up 
to 1.85 GeV with the unitarised chiral expressions. 

In addition, we have also performed chiral fits to the other solutions of Estabrooks et 
al. [^| above 1.5 GeV. It is found that fits with reasonable values for the parameters can 
be obtained for data sets A, B and D whereas we were not able to find an acceptable fit 
for data set C. As far as the x 2 is concerned, the best fit was obtained for data set A 
with x 2 — 1-70, the next best was data set B with x 2 — 1-86 and finally data set D with 
X 2 = 2.57. Therefore, on the basis of these fits, we are unable to confirm that the solution 



found by the group of Aston et al. |39[ is indeed the physical one. Note, however, that in 
ref. |63| arguments were put forward which entail that the solutions A and D of |38 



arc 



unphysical. This result, together with our fits, then favours solution B, in agreement with 
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the findings of ref. |39[]. In what follows, we shall not discuss the additional solutions of 
Estabrooks et al. any further. 

In the region of 1.9 GeV, the experiment by Aston et al. found a second scalar resonance 
in the S-wave channel. If we aim to describe the data up to such energies a second scalar 
nonet should be implemented in the unitarised chiral expressions. To this end we add 
in eqs. Q2.4|) and ( |C. 1| ) to ( |C.5|) a second nonet with degenerate mass M^/ and new scalar 
couplings c' d and c' m analogous to the first nonet. In addition we have also included two more 
vector octets with new couplings Gy> and Gyn. These couplings can be estimated from 
the decay rates K*(U10) -> Kn and if* (1650) -> Ktt @ with the results Gy< = 12 MeV 
and Gyn = 23 MeV, presumably with large errors. For our first fit as in section 4 we fix 
the values of c d and c m to the central values of eq. (|4.3| ). Since we would like to have a 
good description for the 1 = 3/2 channel up to 1.9 GeV, we chose to employ the unitarised 
background expression (|4.7|) as described in section 4. Using the chiral expressions also for 
/ = 3/2, even including the higher scalar and vector resonances, did not provide acceptable 
fits in the region above 1.5 GeV. The fit parameters are then found to be 

M K * = 1.26 GeV, M s , = 2.17 GeV, c' d = 31.6 MeV, c' m = 46.9 MeV, 

c£l = 0.174, c Kv = -0.890, c Kn > = -0.722, (6.4) 

a = -0.760 GeV" 1 , (3 = - 1.300 GeV" 2 , 7 = 3.373 GeV" 4 , 

with a x 2 — 753.3/166 d.o.f. = 4.54. The resulting \ 2 is found to be rather large indicating 
that the estimate of eq. ( |4.3| ) for a and c m does not provide a good description of the data 
if inelastic channels and a second scalar resonance are included. We shall come back to a 
discussion of this point below. Analogous to section 4, the next fit is performed with the 
constraints c m = Cd and c' m = c' d . The resulting fit parameters are 

M K * = 1.29 GeV, M s > = 1.86 GeV, c d = 25.8 MeV , c' d = 14.8 MeV, 

c£l = 0.179, c Kv = 0.504, c Krf = - 1.064, (6.5) 

a = -0.725 GeV" 1 , (3 = -1.817 GeV" 2 , 7 = 4.385 GeV" 4 , 

with a x 2 = 470.9/166 d.o.f. = 2.84. This fit is shown as the short-dashed line in figs. 5 
and 6 for the modulus of the amplitude a and the phase <fio respectively. For comparison, 
in addition, as the dotted line we display the best fit for the elastic case of section 4, 
eq. (TO - 



Our notation for the shown data sets is as follows: solution B of Estabrooks et al. 38 



full circles; solution A of Aston et al. [3~5|, Hfl full squares; solution B of [|3"9"|, f4"0|| empty 
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05° 0-6 




Figure 5: Inelastic fits for clq. The notation for the data sets is as follows: solution B 



of full circles; solution A of [B3,f|0 full squares; solution B of |I§],^0| empty diamonds; 
solution C of |39|,[|(| empty squares; solution D of p9| , |i0| empty triangles. Dotted line: fit 
of eq. (|4. 10Q ; short-dashed line: fit of eq. (|6.5|) ; solid line: fit of eq. ( |6.7|) ; long-dashed line: 
fit of eq. (|6Uq ). 



diamonds; solution C of |39], empty squares; solution D of []39], |4(| empty triangles. 
Although in our fits we have only included data of Estabrooks et al. and Aston et al. up 
to 1.85 GeV where the partial wave analysis of Aston et al. gave a unique solution, in our 
figures all data up to 2.5 GeV are displayed for convenience. Data points where is larger 
than 180° have not been included in the fit since they violate unitarity. This is easily seen 
from the relation 



a sinO ) = ^(l - 7]q /2 cos(25q 



rfj 2 cos(25q 



(6.6) 



which can be derived immediately by taking the imaginary part of eq. (|4.1|) , because for 
< Tjl < 1 the right-hand side is always positive. This shows that solutions B and D 
of [[39|,[40|], also displayed in figs. 5 and 6, are unphysical. 

Finally, we perform a fit leaving both a and c m as free parameters but keeping the 
constraint c' m = c' d because otherwise we would have too many parameters, and the fit 
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Figure 6: Inelastic fits for O - The notation for the data sets is as follows: solution B 
of |38| full circles; solution A of |^,f|0| full squares; solution B of |39],[|0| empty diamonds; 
solution C of |39],fi0| empty squares; solution D of p9|,|40| empty triangles. Dotted line: fit 
of eq. ( |4.10| ); short-dashed line: fit of eq. (|Q|); solid line: fit of eq. (|6.7|); long-dashed line: 



fit of eq. (|O0|) 



would become unstable: 



M KS = 1.30 GeV, M s > = 1.84 GeV, c d = 29.3 MeV, c„ 



8.8 MeV . 



13.1 MeV. 



a 



1/2 

--Kn 



0.151, c Kv = 1.369, c Krf = -1.388, (6.7) 



3.436 GeV" 



- 0.816 GeV" , = -1.045 GeV"' , 7 

with a x 2 = 408.8/165 d.o.f. = 2.48. The resulting amplitude and phase are shown as the 
solid line in figs. 5 and 6. 

As for the fits in the elastic region presented in section 4, also in the inelastic case 
rather different values for the scalar couplings can lead to similar fit quality. Therefore, 
as in section 4, we next try to put further constraints on these couplings to single out 
physically sensible values. As will be discussed in a subsequent publication [pf| , the high- 
energy or short- distance behaviour of the scalar form factor can provide such constraints. 
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Requiring the tree-level form factor in the large- N c limit to vanish at infinite momentum 
transfer, for an arbitrary number of scalar resonances, one finds the following relations: 

4^c mi c di = fl and J^T ( c ™ ~ °di) = °> ( 6 - 8 ) 

i i Si 

where M Sl , c di and c mi are the mass and the couplings of the ith resonance respectively. 

In the case of two resonances, the short-distance constraints of eq. ( |6.8| ) can be used to 
express the couplings of the second resonance c d and c' m in terms of the resonance masses 
and the couplings of the first resonance. For this fit, we find that ck v is zero within the 
fit uncertainties. Thus we decided to set ck v = 0. Fitting the remaining parameters, we 
obtain 

M K * = 1.41 GeV, M s > = 2.69 GeV, c d = 22.8 MeV, 

c m = 26.0 MeV, cjg = 0.178, c Kr f = 0.151, (6.9) 

a = -0.750 GeV -1 , f3 = - 1.553 GeV" 2 , 7 = 4.041 GeV" 4 , 

with a x 2 = 505.6/167 d.o.f. = 3.03. The values for c' d and c' m which resulted from the 
fit of eq. (|6.9p are c' d = 43.8 MeV and c' m = 35.3 MeV. As our last fit, we employ the 
constraints ( |6.8| ) together with the reasonable assumption that c m = Cd and again ck v = 0. 
In this case, the fit parameters are found to be 

M K * = 1.40 GeV, M s > = 2.57 GeV, c d = 23.8 MeV, 

c&l = 0.173, c Kv , = 0.112, (6.10) 

a = -0.758 GeV" 1 , (3 = - 1.473 GeV" 2 , 7 = 3.929 GeV" 4 , 

with a x 2 = 507.1/168 d.o.f. = 3.02. The resulting x 2 is even slightly smaller than for the 
previous fit without the additional constraint c m = c d . For the couplings of the second 
resonance from this fit we obtain c' d = c' m = 39.6 MeV. The latter fit is shown as the long- 
dashed line in figs. 5 and 6. Let us now come to a general discussion about the quality of 
the fits for the inelastic Kn scattering amplitude. 

Even for our best fit of eq. (|6.7| ), the resulting x 2 is larger than two. In our opinion 
the reason for this is twofold. On the one hand, the systematic uncertainties for the data 
above 1.3 GeV seem to be underestimated. This is seen if the data are fit with a Breit- 
Wigner- resonance plus background ansatz [|39| . Even increasing the number of parameters 
in the parametrisation of the background, the lowest x 2 is found to be of the order of 
1.5. On the other hand, additional inelastic channels like Knirn are still missing in our 
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theoretical expressions. This could also lead to deficiencies in the description of the data 
in the higher-energy region. However, in view of our first remark, the experimental data 
are described rather well which in turn implies that additional inelastic channels should 
not be too important. This is supported by the fact that those contributions are of higher 
order in the chiral and 1/N C expansions and should be suppressed. 

Comparing the various fits for the inelastic channel, like in section 4 one finds that 
similar fit quality can be achieved with rather different values for the couplings of the 
scalar resonances. In view of this fact it seemed desirable to put further constraints on 
these couplings from general considerations. These are provided by the short- distance 
constraints ( |6.8| ) which arise if one requires that the tree-level scalar form factor in the 
large-iV c limit should vanish at infinite momentum transfer. For this reason, we believe 
that the fits of eqs. ( |6.5| ) and especially (|6.1(J| ) are physically most sensible, although for 



the description of the data, the fits of eqs. ( |6.7|) and (|6.9|) lead to very similar results. 

As has been mentioned at the beginning of section 5, experimentally it was observed 
that the Kr] channel is not important for S-wave Kn scattering up to 2 GeV. Thus, to 
get a further understanding about which of our fits is most physical, let us investigate 
how sensitive the fits are concerning a removal of the Kr] channel. For simplicity, we just 
consider the fits of eqs. ( |6.5| ) and ( |6.10| ). 



If the Kr] channel is removed in the fit of eq. ( |6.5| ), the x 2 is increased to x 2 = 
664/167 d.o.f. = 3.98 On the other hand, if an analogous fit without the Kr] channel 
and the constraints c m = Cd and c' m = c' d is performed, a fit with similar quality than ( |6.5| ) 
can be obtained: 



M K * = 1.29 GeV, M s , = 1.94 GeV, c d = 27.9 MeV, 

c' d = 20.1 MeV, ell = 0.173, c Ktf = -0.860, (6.11) 

a = -0.846 GeV" 1 , f3 = -0.917 GeV" 2 , 7 = 3.403 GeV" 4 , 

with a x 2 — 479.6/167 d.o.f. = 2.87. Putting back the Kr] channel for these fit parameters 
only increases the x 2 to 522/167 d.o.f . = 3.13. Therefore, as far as the sensitivity with 
respect to the Kr] channel is concerned, the fit (|6.11|) is preferred over the fit 



Performing the same exercise of removing the Kr] channel in the fit of eq. (|6.10j ) , the x 2 
is even slightly decreased to 488/168 d.o.f. = 2.90. In this sense there is some indication 
that the parameter set of eq. ( |6.10| ) is the most physical of our fits to the S-wave Kn 
scattering amplitude. 

In the original publications where the inclusion of resonances in the framework of xPT 
was worked out MM, it was suggested that the order p A coupling constants L 5 and L 8 are 
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saturated by the contribution of the scalar octet resonances to a very good approximation. 
In the case of two resonances these contributions are given by: 



CdCm c d C m 



„/2 



Ml 



Ml 



and 



2M| 



+ 



2M|, 



(6.12) 



Thus, assuming resonance saturation of L 5 and L 8 and using our fit parameters we are 
in a position to calculate these chiral couplings. Taking the fit of eq. ( |6.10| ) including the 
short-distance constraint as an example, we then find 



0.52 • 10" 



and 



0.26-10" 



(6.13) 



in good agreement to the most recent determination of these constants |46],|47j and also 
supporting lower values than were found in earlier estimates. The fact that in this approach 
we obtain L5 = 2L§ follows immediately from the second constraint of eq. ( |6.8| ). Similar 
values for L 5 and L 8 are also obtained for the other fits to the inelastic scattering amplitude. 



7 The T-matrix poles 

In the energy region above the Kir threshold and below roughly 2 GeV, the possible 
presence of three scalar resonances in the S-wave I — 1/2 channel has been discussed 
in the literature. The only one which is clearly observed in the experimental data is the 
i^o(1430) P^ , |55| . Less obvious is the existence of the i^o(1950) resonance JHJ in the region 
around 2 GeV, and concerning the third, the lowest-lying k meson [p!9| , |20| , [6^ -|69l , there has 
been some controversy in recent years |)4],[7(|[7T] . 



We now address this problem by searching for possible poles of the matrix £(s), eq. 
associated with these resonances in the nearby unphysical Riemann sheets. As discussed 
above, the influence of the Kt] threshold in the S-wave 1 = 1/2 partial wave amplitude 
was found to be very small and hence, in order to simplify the notation, we shall label 
the different Riemann sheets as if we had only a two channel problem with and Krj'. 
In this way, the first or physical Riemann sheet is indicated by (+, +), the second one by 
(—,+), the third one by (— , — ) and the fourth by (+,—). The signs between brackets 
refer to the sign ambiguity present in the definition of the modulus of the centre-of-mass 
three-momentum due to the presence of the square root of eq. (|3.5| ). In this way, the '+' 
sign indicates that the imaginary part of the modulus of the three- momentum is always 
positive in the whole s-plane and the '— ' sign refers to the contrary. The sign between 
brackets on the left side corresponds to the Kir system and the one on the right side to the 
Kr]' state. We start by investigating the pole positions for our preferred fit of eq. ( |6.1U| ). 
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At the end of this section, we shall also compare with the pole positions found for some of 
the other fits obtained in sections 4 and 6. 

The nearby Riemann sheet for the k resonance is the second one, since this state should 
appear with a mass much lower than the threshold of the Krj' channel. In fact, we find a 
pole in the second Riemann sheet (— , +) at 

^ k = (708 - % 305) MeV . (7.1) 

This pole position is very similar to the one found in ref. |21J at (779 — i 330) MeV.Q Note 



that the main differences between ref. [20|, applied to the specific case of S-wave 1=1/2 
scattering, and the present work are that in the former unphysical cut contributions were 
not included, the Krf channel was not considered (although its influence below 1.3 GeV 
is very smooth), and the short distance constraints were not imposed. However, in the 
present work, although we have included all this information, the pole ( |7.1|) corresponding 
to the k resonance appears in a similar position to the one found in ref. [J20| . (We are 
talking about a resonance with a width of about 600 MeV; hence small variations in its 
mass and width under realistic changes in the dynamical input are completely natural.) 



These findings are expected since in ref. P0[ it was estimated that the contribution from 



the unphysical cuts in the physical region below 0.8 GeV amounts just to a few percent. 
Thus the stability of the expansion in the physical region in terms of the unphysical cut 
contributions at these energies seems to be justified. 

We have checked that this pole is of dynamical origin, that is, it is generated through 
the strong rescattering (unitarity loops) of the Kir channel. For instance, this pole survives 
when one removes all the explicit tree-level resonant contributions, as was already estab- 
lished in ref. [pOR . Because of this fact and also because the n is a very wide resonance, 



which implies that its contribution to the physical data is rather soft, one can expect that 
the existence and properties of this pole are rather model dependent. This is basically 



what is observed in the recent work ||71|| , where a different conclusion has been obtained. 
However, contrary to what is concluded in this reference, the model dependence does not 
imply that the k does not exist. It means that a good dynamical model, based on first 
principles, is required, in order to assess the question whether there is or is not a pole to 
be associated with the k resonance. This is in fact one of the aims of the present work. 

Before finishing the general discussion with respect to the k, let us remind that this 
pole, together with the <r(500), a (980) and a strong contribution to the / (980), gives 

4 Employing the pole position one can identify the real part with the mass of the resonance and the 
modulus of the imaginary part with one half of its width. For narrow objects, this identification agrees 
with the Breit-Wigner resonance picture. 
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rise to the lightest scalar nonet ||20|| . In this work, the whole J = ++ nonet was found 



to be of dynamical origin, that is, due to the strong rescattering in this channel between 
the pseudoscalars from the lowest order amplitudes. In fact, taking the results of 
ref. [0, it is possible to go to the SU(3) limit by considering the limit of equal masses for 
the pseudoscalars. In that way, one can see how the poles move continuously giving rise 
to an octet of degenerate scalar resonances plus a singlet. 

Now we turn to the -K"q(1430) resonance. This state is clearly seen as a steep rise in 
the phase, figs. 3 and 6, and as a peak structure in the modulus of the amplitude, figs. 2 
and 5. The threshold of the Krj' is very close to the expected mass of this resonance and 
hence we should investigate the pole positions both in the second (— , +) and in the third 
(— , — ) Riemann sheet, which are the nearby ones in this case. We find the following pole 
positions for the fit ( |6.10 ): 



Sheet Pole Position 

(- +) v^ . ( i430) = ( 1450 - i 142 ) MeV ( 7 - 2 ) 
(" -) v^uso) = (1358 - i 243) MeV 

There is a clear dependence in the position of the pole with respect to the sheet. This is 
expected since the resonance, although originating from a tree level pole included around 
1.4 GeV, has strong unitarity corrections and is also influenced by the nearby Krj' threshold. 
In any case, the pole position found in the second Riemann sheet is very close to the effective 
values of the mass 1429 ± 4 ± 5 MeV and width 287 ± 10 ± 21 MeV given in refs. |9 



deduced from a simple Breit-Wigner fit to the experimental data. It is worth to stress that 
this is the lightest preexisting 1 = 1/2 scalar state. 

Finally, some structure in the phase and modulus of the amplitude around 1.7—1.9 GeV 
is also observed in the experimental data by Aston et al. (solutions A and C of ref. 
in figs. 5 and 6, reminding that solutions B and D of this reference are unphysical). In 



ref. [39], an additional simple Breit-Wigner plus background fit was performed for this 
energy region only, establishing the presence of the (1950) resonance. The effective mass 
and width given in that work varied as a function of the chosen solution for the partial wave 
amplitude. For solution A, a mass of 1934 ± 8 ± 20 MeV and a width of 174 ± 19 ± 79 MeV 
was found, whereas for the unphysical solution B the mass was 1955 ± 10 ± 8 MeV and the 
width 228 ± 34 ± 22 MeV. As explained above, our fits were performed including data from 



ref. |39] up to 1.85 GeV. In this energy interval the solution of ref. [^] is unambiguous 



and in fact our fits also give rise to a pole in the energy region around 1.9 GeV. For the fit 
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( |S.10| ), the poles in the nearby third (— , — ) and in the second Riemann sheet (— , +) are: 

Sheet Pole Position 



^(1950) = (1731 -i 147) MeV (7.3) 
(-+) ^(1950) = ( 1908 - ^27)MeV 

One can see a large variation in the position of this pole from one sheet to another indicating 
a strong departure from the simple Breit-Wigner resonance behaviour. Furthermore, we 
have checked that this resonance is an interplay between the underlying dynamics discussed 
in the previous sections and rescattering effects.^ On the other hand, while the mass from 



the pole position in the second sheet is closer to the effective one of ref. [p9], the width 
from the pole in the third sheet is closer to the one given in that work. 

Let us now come to a comparison of the pole positions resulting from different fits of 
sections 4 and 6. In table 1, we give the poles corresponding to the previously discussed 



resonances in the second sheet for the fits: ( |640D , flOlD , Q, ([L4|), O and fl4lCp . 



Since the last three fits correspond to the elastic case of section 4, only the resonances k 
and _K"q(1430) are considered for them. In this table, we also show the residua of the poles 
associated with the different channels. They are defined by the limit: 



n = lim JKs-s^Tul (7.4) 

with i referring to the considered channel (i = 1,2,3 <-> Kit, Kr], Kr]') and Sr is the 
location of the pole for the resonance in question. The residua presented in table 1 show 
that the k couples strongly to the Kw state and much less to the Kr] and Krf channels. 
The _K"q(1430) and _K"q(1950) couple strongly with the Kir and Kr]' states with comparable 
strengths in both cases. On the other hand, the coupling of these resonances with the Kr] 
state is much weaker. 

When comparing the different fits in table 1, one sees that the properties of the k meson 
are rather stable when changing from one fit to another. For the _K"q(1430) one observes a 
common trend in the inelastic fits, as well as in the elastic one of eq. ( f4.10| ), to give a pole 
around (1.44 — i 0.15) GeV. There is, however, a large difference between the pole position 
and the residua r 1 corresponding to the elastic fit ( [4.9|) and the ones given by the inelastic 
fits. On the other hand, the values from the fit ( |4.4j ) lie in an intermediate region. We will 
further discuss these points after introducing table 2 below. 

For the _K"q(1950) the situation is more unstable when comparing different fits, particu- 
larly with respect to the imaginary part of the pole position in the second sheet. However, 
5 The bare parameters of the second resonance change considerably for the different fits and in fact for 



the fit (3.10) the bare mass is as high as 2.57 GeV. 
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Sheet :(-,+) 


K 


lift (1430) 


Ko(1950) 


(MOD 

r l,2,3 


(708 - i 305) MeV 
4.51, 2.36, 1.68 


(1450 -z 142) MeV 
4.86, 1.07, 4.12 


(1908 -i 27) MeV 
5.01, 0.73, 5.97 


(6.11 

r l,2,3 


) 


(694 - % 329) MeV 
4.83, 2.78, 1.72 


(1447-il60) MeV 
5.32, 1.29, 4.71 


(1966 -i 50) MeV 
5.75, 0.99, 7.03 


G3) 

r l,2,3 


(714-i301) MeV 
4 31 1 55 1 19 

4:.iJl 5 A..OO^ A..A.C/ 


(1429-il46) MeV 
4 8? 85 9 67 


(1880 - i 112) MeV 
4 76 91 4 90 


o 

n 


(701-i320) MeV 
4.79 


(1437 - i 70) MeV 
3.53 


Elastic Case 


H) 

n 


(709 - i 306) MeV 
4.74 


(1386 -i 48) MeV 
2.85 


Elastic Case 


(|4.10 

ri 


) 


(684 - % 279) MeV 
4.10 


(1407 - i 140) MeV 
4.21 


Elastic Case 



Table 1: A collection of the pole positions and residua for the fits of eqs. (|6.10| ), Q6.ll]) , 
(|6.7D, (|4.4|), (|4.9|) and (|4.10|) in the second Riemann sheet. The r» are given in units of M n . 



as can be seen in figs. 5 and 6, all the inelastic fits result in rather similar curves. This 
clearly implies that the role of the ^(1950) is hidden by a very large background and 
hence one needs a very precise approach for such large energies around 2 GeV in order 
to extract the resonance-pole parameters. This is clearly out of the scope of the present 
paper because for such high energies, additional multiparticle states can play a significant 
role and also the contributions of the unphysical cuts become increasingly large. Let us 
note that for the fit ( |6.7|) the pole position corresponding to this resonance is very similar 
to the mass and width obtained in ref. |[39|| , the only reference considered in the Review of 
Particle Physics @ for the #5(1950). 

Table 2 is analogous to table 1, but now considering the third Riemann sheet (—,—). 
Therefore, in this second table, we do not show the k resonance and we only consider 
inelastic fits. The third sheet is in principle the closest one to the physical s-axis when one 
is well above the Kr{ threshold. This is clearly the situation for the 1Tq(1950). For the 
i^o(1430) one also needs to consider the second sheet (— , +) discussed above, since we are 
very close to the Kr( threshold. From table 2, one infers that for this latter case the pole 
associated with the 1^(1430) in the third sheet is very broad, but with similar residua to 
the one in the second sheet, and thus one would expect that the pole driving the effects 
of the ^(1430) on the physical axis should correspond to the one in the second sheet. In 
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Sheet:(— , — ) 


if *(1430) 


fsT *(1950) 


d6.ic|) 


(1358 -i 243) MeV 


(1731 - i 147) MeV 




5.56, 1.65, 4.41 


3.67, 0.24, 4.84 


(6.11) 


(1299 -i 255) MeV 


(1700 -i 151) MeV 


r l,2,3 


5.52, 1.44, 4.07 


3.79, 0.24, 4.62 


(Q 


(1394 - i 188) MeV 


(1772 - i 172) MeV 


r l,2,3 


5.18, 0.84, 3.19 


4.24, 0.58, 4.42 



Table 2: A collection of the pole positions and residua for the inelastic fits of eqs. ( 6.1C ), 
( |6.11| ) and (|6.7|) in the third Riemann sheet. The are given in units of M v . 



fact, the masses and widths of the -K"q(1430) derived from the second sheet poles listed in 
table 1, corresponding to the inelastic fits and also to the elastic one ( |4.10| ), are in fairly 
good agreement with the effective values of ref . |39| . 

Concerning the ^(1950) resonance poles, there is no clear and simple criterion in order 
to decide which pole, the one in the second or in the third sheet, is the leading one. In fact, 
when comparing with the experimental data of solutions A and C of ref. ||39|| , fig. 5, one 
observes a dip in the modulus of the amplitude around 1.7 GeV. This precisely corresponds 
to the real part of the pole position corresponding to the Kq(1950) in the third Riemann 
sheet, table 2. One the other hand, when looking at the phase shifts for the same solutions 
in fig. 6, one sees a local maximum around 1.9 GeV, which coincides with the real part of 
the poles in the second sheet, table 1. Furthermore, as discussed above, in the region of 
the ^(1950) there is a large background that could mimic the effects of the broad and 
lighter pole found in the third sheet. 



8 Conclusions 

S-wave scattering of the pseudoscalar Goldstone bosons is a good testing ground for our 
understanding of strong interactions in the non-perturbative region. At low energies, the 
Goldstone dynamics is constrained by the chiral symmetry properties of QCD. The 
Lagrangian contains all relevant information on the infrared singularities (thresholds, cuts, 
etc.), which allows to reconstruct the scattering S-matrix elements up to local subtraction 
polynomials. The structure of these local terms is also governed by chiral symmetry, but 
their coefficients, the chiral couplings, encode the short- distance information. 

The xPT predictions can be extrapolated to higher energies by taking into account 
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the leading massive singularities, i.e. incorporating the lightest resonance poles through 
a chiral-symmetric Effective Field Theory with resonance fields as explicit degrees of free- 
dom RPJ. Combined with large- N c arguments to organise a well-defined perturbative ap- 
proach, this procedure provides a good understanding of many low and intermediate energy 
phenomena. However, it is not enough to achieve a correct description of the scalar sector, 
because the S-wave rescattering of two pseudoscalars is very strong, making necessary to 



perform a resummation of chiral-loop corrections in order to satisfy unitarity [p!7| - |21| . 

In this paper, we have tried to describe the S-wave Kit scattering up to 2 GeV, through 
a unitarisation of the resonance chiral Lagrangian predictions. In the elastic region a 
reasonably good description can be achieved for the dominant J = 1/2 amplitude, which 
satisfies all low-energy and high-energy constraints. However, the more exotic 1 = 3/2 
sector remains problematic, in spite of the clear improvement provided by the unitarisation 
procedure. Being pure background (no resonances), and with very large contributions 
from crossed channel dynamics, it is more difficult to get a good high-energy extrapolation 
for the 1 = 3/2 amplitude since many missing small corrections could be numerically 
important. Nevertheless, the present discrepancies among different sets of data prevent us 
from reaching any firm conclusion; in fact, the older 1 = 3/2 data appears to be better 
described within the chiral framework. 

Above roughly 1.3 GeV, the 7 = 1/2 amplitude contains inelastic contributions. We 
have incorporated the leading two-body Krj and Krf modes, through a coupled channel 
analysis. The necessary input are the five independent Tk-k^k^ and Tk^^k^ amplitudes, 
which we have computed at the required order. 

The inelastic fit faces the problem of a very unsatisfactory experimental status. The 
data only determines the absolute value of the scattering amplitude, leaving 2 L discrete 
ambiguous solutions to the partial wave decomposition with L being the maximal angular 
momentum contributing at a given energy. Moreover, different experiments give rise to 
different solutions making rather unclear the comparison with theory. Even the masses 
and widths of the dominant scalar resonances are not clearly established. 

We have combined the known theoretical ingredients in order to see which experimental 
solutions are preferred. Up to 1.85 GeV, the Aston data [FMBOJ corresponds to solution B 



of Estabrooks |j38[ . Arguing that this should be the physical solution, we have performed 
several fits in the respective region and studied their behaviour at higher energies where the 
data ambiguities are more severe. Although some of the solutions of ref. |3£§ for energies 
higher than 1.85 GeV are shown to be unphysical, since they violate unitarity and can be 
discarded, the large experimental errors do not allow to unambiguously select the physical 
one. 
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A similar fit quality could be achieved with rather different values of the input param- 
eters. Nevertheless, when imposing the high-energy constraints discussed above, definite 
values for some of the parameters can be given. In addition, assuming saturation of the 
low-energy constants L5 and L$ by scalar-resonance contributions, we were able to estimate 
values for these parameters which are compatible with the very recent work |46|,^7| . We 
have also used our fits to study the position of the T-matrix poles in the complex plane, 
aiming to determine the mass and width of the dominant scalar resonances. 

Clearly, better experimental data are needed. Our results provide a convenient theoreti- 
cal framework to be used in future experiments for analysing the data and resolving possible 
ambiguities. Even with all present shortcomings, our analysis considerably improves the 
knowledge of the S-wave Kir scattering amplitude. In a forthcoming publication |28 



we 



will use this information to perform a detailed investigation of the Kit scalar form factor 
up to 2 GeV. We expect that this could help to improve the determination of the strange 
quark mass from QCD sum rules for the strangeness-changing scalar current p2|-p^1. 
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Appendices 



A Loop functions 

For completeness, we tabulate here the one-loop functions || appearing in the unitarity 
chiral correction of eq. (2.5) . The basic integral is given by: 



■Jpq(s) 



IQtt 2 Jo 



1 



32tt 2 



2 + 



dx log 
'A 



Ml - sx{l - x) - A PQ x 



Ml - A PQ x 



(A.l) 



PQ 



^PQ, 



log 



Ml, 



log 



(s + Xpq) 2 - A 



PQ 



fs — A 



PQ) 



A 



PQ 



with 



Sp Q = Ml + M z Q - A P q = M l P ~M l Q - \ P q 
From Jpq(s), one defines the related quantities: 

Jpq(s) = Jpq(s) - 2kp Q {(i) ; K PQ (s) 



s _ ( Mp + M Q ) 2 ] \s - (M P - M Q f 



ArQ J PQ (s) ; L PQ (s) = J PQ (s) ; (A.2) 



and 



PQ 



M PQ (s) = — (s-2Ep Q )Jp Q (s)+ , w 



2s 



Jpq(s) - Sj'f 



PQ 



Is 



- l kpM + 28^ 



(A.3) 



Here, 



327T 2 Ap Q 
1 



Ml log (M|//i 2 ) - M 2 Q log (M 2 Q //r 



E PQ , n MlMr 



V = - t ~ J v 4- 2- 



A 3 
^PQ 



log 



Mp, 



(A.4) 



and pi denotes the renormalisation scale. 
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B 77-77' mixing 

Let us consider the two-dimensional space of isoscalar pseudoscalar mesons. We collect 
the SU(3) octet and singlet fields in the doublet 77^ = (778,771). The quadratic term in the 
Lagrangian takes the form 

C = ^d^llCd^B-^M'vB, (B.l) 



with 



5 1 + ei / I M 2 8 M 2 



The coupling to the scalar resonances generates kinetic matrix, 

8c d c m 02 8c d c m 02 2v /2 8c d c m 

68 = 7^m| ' 61 = 7^ Ml ' 5 = - — jHf! AK - (R3) 

and modifies the lowest-order values of the mass-matrix elements: 
Ml = Ms +^l( M ' + l M K^), 

Ml = M%+ Ml +^{ M t + \Ml^y (B.4) 



2^2 / 8c 2 



Ml = 7T~ Ak-it ( 1 + -pTTTo %Mf 



PM S 

i 

18 " ^Kn \^J2M2^K 

02 

Here, Mq denotes the U(1)a anomaly contribution to the 771 mass, Mj (i = 1,8) the 
0{p 2 ) contributions to the singlet and octet isoscalar masses, 

02 1/02 2 \ 02 1/02 2 \ 

M 8 = -(4 M K M w J, M 1 = -[2 M K + M„ J, (B.5) 

02 02 02 02 

and Ax7r= M K — M n , with M n and the pion and kaon masses at 0(p 2 ) in the 

chiral expansion. 

To first order in ei, es and 5, the kinetic matrix /C can be diagonalised through the 
following field redefinition: 



Vb 



VS \ r7\l2 T V" rvl/2 T ryl/2 ( ^ ~ ^ ~ 



Z l/2 . f, = Z l/2 . '/« 1 Z V2 J . £ . Z l/2 = j Z l/2 



2 2 



* r v -f !-t 



(B.6) 
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In the fj basis the mass matrix takes the form 

M = Z l ' 2T ■ M ■ Z 1 ' 2 , (B.7) 

where 

Mi ) + + 8e " (c ^ - Cj) + 5 MjA K -,) , (B.8) 

S " s 

( i i 8c m (c m — Q) 2 ^ C mCd A/ r2\ 

K V + PMI ZMk pmI m °) ■ 

The physical mass eigenstates are obtained by diagonalising the matrix Xi with an 
orthogonal transformation, 

T-i ^ „ 7^ t-> ( V \ ^ I cos 6 P — sin6>p \ 

\ V ) \ sm&p cos Op I 

(B.9) 

One gets the relations 

M 2 , + M 2 = M? + M£, Ml -Ml = \J (Mi — M|) 2 + 4Mi 8 , (B.10) 



M 8 2 


2 

= M 8 +■ 


M 2 


= M 2 (l- 


Ml8 


2v/2 
3 z 



and 



tan, P = = M *>7^ = ^ = -3-. (B.ll) 

M 2 8 M 2 8 M 2 -Mf Ml -M^ 
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C Inelastic channels 



Rr 2 A 2 



1 [ 2MI + MI 2c m c d A K ^ 
V2p K \ 3 + 3fMl {M P ^ + 1K) 3f 2 M 2 

C d S + (c m - C d )H Kn 



+ 



p 



+ 



+ 



Mi* - s 
c d u + (c m - c d )T, Kv 



Mi* - 



c d t + 2(c ri 



c d )M 2 K 



Ml 



c d (s - V \ - M 2 K ) + 2c m M 2 K 
c d (u-p 2 v - M 2 K ) + 2c m M 2 K 
c d (t-p 2 -M 2 ) + 2c m M 2 



(C.l) 



_L J Q/ _ „2 _ § 2 1 2 , 4c m Q Afly 2 2 

4/|.\^ P " 3 Mk 3 IK+ 3PM 2 ^ t + 2p r> IM «> 3 PMI 



+ 



c d s+ (c m - C d )^ Kv 



+ 



+ 2 



Mi* - s 
c d u + (c m - c d )£ K7r 



Mi* - u 
+ 2(c m - c d )M 2 K 



M„ 2 -t 



c m (3M W 2 - 5M|) - c d (s - p 2 ~ Ml) 
c m (3M 2 W - 5M 2 K ) - Cd (u-p 2 v - M 2 K ) 



2c m M 2 + Cd (t-p 2 - Ml) 



+ 



3G 2 , \ s(t -u)- (p 2 - M 2 K )A Kn u(t -s)- (p 2 - M 2 K )A Kn 



P 



Ml* - s 



Ml* - u 



(C.2) 



Tk^i^Ktu — 1 2 1V±K 



I M 2 Ac mC d A 2 Rn 



+ 



P 



+ 



+ 



M 2 - s 



M 2 K * - u 

1 



+ 



2 c m A^ 

3 Ml -t 



3f 2 M 2 

c d (s-k 2 ,- M 2 K ) + 2c m M 2 K ] [c d (s - k' 2 - M 2 K ) + 2c m M 2 K 
~c d (u-k 2 - M 2 K ) + 2c m M 2 K ] [c d (u - k' 2 - M 2 K ) + 2c m M 2 K 
c d t + 2(c m - c d )M 2 K ] [c d (t-k 2 - k' 2 ) + | c m (2M 2 K + M 2 
c d t + 2(c m -c d )M 2 K ]^ j, (C.3) 
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T, 



V2_ 
3/i 



Ml - \Ml 



C m Cd L^Ktx 



4dA? 



|^(9t - 3k 2 + 3^ + 16A K7r ) + 

S J 



m Kit 



M 2 



+ 



P 



+ 



+ 



M|-* - s 



M 2 s -t 



+ 



16 c m A^tt r 
3 M 2 Si -t 



c d (s-k 2 - M 2 K ) + 2c m M 2 K ] [c d (s - k' 2 - M 2 K ) + c m (5Ml - 3M 2 ) 

c d (u - k 2 - Ml) + 2c m M 2 K ] [c d (it - k' 2 - M 2 K ) + c m (5M 2 K - 3M 2 

c d t + 2(c m - c d )M 2 K ] [c d (t-k 2 ,- k' 2 ) + | c m (4M 2 K - M 2 
c d t + 2(c m -c d )Ml\]\ , (C.4) 



T; 



Kr/$—>Kr]$ 



m 12 



3f 2 M 2 



-(M + 32A Kn ) + 



(9t - 3k 2 - 3k' 2 - 2M 2 ) • 2CmCd AKn , oo A N , 14 <4 A L 



f 2 M 2 s 



1 



f 2 lM 2 -s 



M 2 ,-u 



c d {s-k 2 v - M 2 K ) + c m {5M 2 K - 3M 2 )} [c d (s - k' 2 - M 2 K ) + c m {5M 2 K - 3M 2 ) 



c d (u-k 2 - Ml) + c m (5M 2 K - 3M 2 ) c d (u - k' 2 - M 2 K ) + c m (5M£ - 3M 2 



+ 



M 2 s -t 



c d t + 2(c m - c d )M 2 K c d (t-k 2 - k' 2 ) + | c m (8M 2 K - 5M 2 
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+ M2 ~ — [c d t + 2{c m - c d )M 2 K \ [c d (t-k 2 ,- k' 2 ) + f c m (AM 2 K - M 2 ) 
9 Gy \s(t -u) + (k 2 - M 2 K )(k' 2 - M 2 K ) + u(t -s) + (k 2 - M 2 K )(k' 2 - M 2 K 



+ 



2f 



Ml*-s 



Ml, - u 



(C.5) 



Several remarks concerning the above expressions are in order. p v , k v and k' v are 
momenta of the rj^s. In eqs. ( |C3| ) to ( |U.5| ), k v refers to the initial state r\i and k!^ to the 
final state rjj. For example in the calculation of the amplitude Tk v ^kti', we have k 2 = M 2 
and k' 2 = M 2 ,. 
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D Input parameters 

For the convenience of the reader, in this appendix we collect the values of all constant 
input parameters used in our numerical analysis: 

M w = 138.0 MeV, M K = 495.7 MeV, 

M v = 547.3 MeV, M„, = 957.8 MeV, 

f n = 92.4 MeV, f K = 112.7 MeV, 

M p = 770.0 MeV, M K * = 896.1 MeV, 
M v , = 1440 MeV, M v „ = 1710 MeV, 

G v = = 65.3 MeV, 
G v , = 12 MeV, Gyn = 23 MeV, 
M s = M Sl , M s , = M s[ , 
sin9 = -1/3 w 20°. 
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